COSMOS: 

3D weak lensing and the growth of structure 
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ABSTRACT 

We present a three dimensional cosmic shear analysis of the Hubble Space Telescope COSMOS 
survey*, the largest ever optical imaging program performed in space. We have measured the 
shapes of galaxies for the tell-tale distortions caused by weak gravitational lensing, and traced 
the growth of that signal as a function of redshift. Using both 2D and 3D analyses, we measure 
cosmological parameters f2 m , the density of matter in the universe, and as, the normalization 
of the matter power spectrum. The introduction of redshift information tightens the constraints 
by a factor of three, and also reduces the relative sampling (or "cosmic") variance compared 
to recent surveys that may be larger but are only two dimensional. From the 3D analysis, 
we find as (f2 m /0.3) ' 44 = 0.866^Q;ggg at 68% confidence limits, including both statistical and 
potential systematic sources of error in the total budget. Indeed, the absolute calibration of 
shear measurement methods is now the dominant source of uncertainty. Assuming instead a 
baseline cosmology to fix the geometry of the universe, we have measured the growth of structure 
on both linear and non-linear physical scales. Our results thus demonstrate a proof of concept 
for tomographic analysis techniques that have been proposed for future weak lensing surveys by 
a dedicated wide-field telescope in space. 
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1. Introduction 

The observed shapes of distant galaxies be- 
come slightly distorted as light from them passes 
through foreground mass structures. Such "cos- 
mic shear" is induced by the (differential) grav- 
itational deflection of a light bundle, and hap- 
pens regardless of the nature and state of the 
foreground mass. It is therefore a uniquely pow- 
erful probe of the dark matter distribution, di- 
rectly and simply linked to theories of struc- 
ture formation that may be ill-equipped to pre- 
dict the distribution of light ( for reviews, see 
Bartelmann fc Schneider! 120001 : IWittmanl 120021 
Refregierll2004l) . Furthermore, the main difficul- 
ties in this technique lie within the optics of a 
telescope that has been built on Earth and can 
be thoroughly tested. It is not limited by system- 
atic biases from unknown physic s like astrophys- 
ical b i as (iDekel fc Lahav 1999 ;_ Hoekstra et al 



shear measurement methods iBridle et all 12001 



Rcfregier fc Baconl[ 2003; Bernstein & Jarvra l200 
Massev fc Refregierl200alMandelbaum et al.ll200 



12: 

05; 



Kuiikeij2006tlNakaiima fc : Bernsteiij2006tlMassev et al 



2006). With several ambitious plans for dedi- 
cated telescopes both on the ground (e.g. CTIO- 
DES, Pan-STARRS, VISTA/ VST-KIDS, LSST) 
and in space (e.g. DUNE, SNAP, and other possi- 
ble JDEM incarnations), the importance of weak 
lensing in future cosmological and astrophysical 
contexts seems assured. 

In this paper, we present statistical results 
from the first space-based survey comparable 
to those from dedicated ground-based observa- 
tions. The "cosmic e volution survey" (COSMOS; 



2002ri ISmith et all 120031 : IWeinberg et all l2000h 



lected galaxy clusters ( 


Huterer & White! 20031: 


Pieroaoli. Scott & White 


20011 Viana. Nichol & Liddle 


20021). 



The study of cosmic shear has rapidly pro- 
gressed since the simultaneous detection of a co- 



Scoville et al.l [2007a)) combines the largest con- 
tiguous expanse of deep imaging from space, 
with extensive, multicolor follow-up from the 
ground. High resolution imaging is particularly 
needed for weak lensing because the shapes of 
galaxies that would also be detected from the 
ground are much less affected by the telescope's 
PSF, and a much higher density of new galaxy 
shapes are resolved. This allows the signal to 
be measured on smaller physical scales for the 
first time. Parameter constraints from our sur- 



vey still carry a fair deal of statistical uncertainty 

herent signal by four independent groups ([Bacon, Refregier^^Us cosmic variance in the finite survey size; 



Van Waerbeke et al 



20001 iKaiser et al.l l2000l: IWittman et all 1200a 



2000). Large, dedicated 



surveys with ground-based telescopes have re- 
cently measured the projected 2D power spec- 
trum of the large-scale mass distribution and 
drawn competit ive constraints o n cosmologi- 
cal p a rameters (Brown et al. 20031 Bacon et al. 
20031: lHamana et all 120031: IJarvis et all 1200a 



based surveys (Rhodes. Re 


Eregier & Grothl l200lt 


Refregier Rhodes. & Groth 


20021 Rhodes et all 


2004: Hevmans et al.l 2005) 


More importantly, 



Van Waerbeke. Mellier fc Hoekstral2005l:lMassev et al 



20051 : iHoekstra et al.l 120061 ). The addition of 
photometric redshift estimation for large num- 
bers of galaxies has led to the first measure- 
ments of a ch anging lensing signal as a func- 
tion of redshift dBacon et al.ll2005b lwittman|[2005l 



Semboloni et alj|2006h 

The shear measurement methods used for 
these ground-based surveys have been precisely 
calibrated on simulated images containing a 
known shear signal by the Shear TEsting Pro- 
gram (STEP; iHevmans et"aH 120061 ; iMassev et al 



2007h . This program has also sped the develop- 
ment of a next generation of even more accurate 



the potential level of observational systematics 
is much lower from space than from the ground, 
where the presence of the atmosphere fundamen- 
tally limits all weak lensing measurements. 

Extensive ground-based follow-up in multiple 
filters has also provided photometric redshift es- 
timates for each galaxy. Lensing requires a purely 
geometric measurement, so knowledge of the dis- 
tances in a lens system as well as the angles 
through which light has been deflected are es- 
sential. We have extended cosmic shear anal- 
ysis into the information-rich three dimensional 
shear field. Our constraints on cosmological pa- 
rameters are tightened by observing independent 
galaxies at multiple redshifts, and the separate 
volume in each redshift slice reduces the cosmic 
variance. Furthermore, we can directly trace the 
growth of large-scale structure on both linear and 



2 



non-linear physical scales. Although these re- 
sults are still limited by the finite size of the 
COSMOS survey, they provide a "proof of con- 
cep t" for tomographic techniques suggested (by 
e.g. Taylor 20021: iBernstein k. J ami 12004 : iHeavensl 
2005t iTavlor et al.l 120061 ) for future missions ded- 
icated to weak lensing. Throughout this paper, 
we have assumed a flat universe, with the Hubble 
parameter h = 0.7. 

This paper is organized as follows. In we de- 
scribe the data and analysis techniques. In we 
present a traditional 2D "cosmic shear" analysis of 
the two-point correlation functions, demonstrat- 
ing the level to which systematic effects have been 
eliminated from the COSMOS data. In 21 we ex- 
tend the analysis into three dimensions via redshift 
tomography. We show how the signal grows as a 
function of redshift, and directly trace the growth 
of structure over cosmic time, on a range of phys- 
ical scales. In EJSJ we use the measured statistics 
from both the 2D and 3D analyses to derive con- 
straints on cosmological parameters. We conclude 



2. Data analysis methods 
2.1. Image acquisition 

The COSMOS field is a contiguous square, 
covering 1.64 square degrees an d centered at 



10:00: 2 8.6, +02:12:21.0 (J200 0) (|Scoville et al 
2007bl : iKoekemoer et al.l 120071 ). Between Octo- 



ber 2003 and June 2005, the region was com- 
pletely tiled by 575 slightly overlapping pointings 
of the Advanced Camera for Surveys (ACS) Wide 
Field Camera (WFC) with the F81AW (approx- 
imately /-band) filter. Four, slightly dithered, 
507 second exposures were taken at each pointing. 
Compact objects can be detected on the stacked 
images in a 0.15" di ameter aperture at 5 <7 down 



to F8UW AB = 26.6 (|Scoville et al.ll2007al) 



The individual images were reduced using the 
standard STScI ACS pipeline, and combined 



using MultiDrizzle (jFruchter fc Hookl [2002; 



Koekemoer et al.l 120071 ) . We took care to opti- 



mize various MultiDrizzle parameters for precise 
galaxy shape measu rement in the stacked images 



(|Rhodes et al.l 120071) . We use a finer pixel scale 



of 0.03" for the stacked images. Pixelization acts 
as a convolution followed by a resampling and, 
although current algorithms can successfully cor- 



rect for convolution, the formalism to properly 
treat resampling is still under development for 
the next generation of methods. We use a Gaus- 
sian DRIZZLE kernel that is isotropic and, with 
pixfrac= 0.8, small enough to avoid smearing 
the object unnecessarily while large enough to 
guarantee that the convolution dominates the re- 
sampling. This process is then properly corrected 
by existing shear measurement methods. 

2.2. Shear measurement 

The detection of objects and measurement of 



their s hapes is fully described in lLeauthaud et al 



(200 7|). Mode l ling o f the ACS PSF is discussed in 
Rho des et al. I (l2007h . Here we provide only a brief 
summary of the important results. 

Objects were detected in the reduced ACS im- 
ages using SourceExtractor ( Bertin fc Arnoutsl 
1996). To avoid biasing our result, the detection 
threshold was set intentionally low: far beneath 
the final thresholds that we adopt. The catalog 
was finally separated into stars and galaxies by 
noting their positions on the magnitude vs peak 
surface brightness plane. Objects near bright stars 
or any saturated pixels were masked using an au- 
tomatic algorithm, to avoid shape biases due to 
any background gradient. The images were then 
all visually inspected, to mask other defects by 
hand (including ghosting, reflected light and as- 
teroid/satellite trails). 

The size and the ellipticity of the ACS PSF 
varies over time, due to the thermal "breathing" 
of the spacecraft. The long period of time dur- 
ing which the COSMOS data was collected forces 
us to consider this effect. Although other strate- 
gies have been demonstrated successfully for ob- 
servations conducted on a shorter time span, it 
would be inappropria te for us to assume, like 
Lombardi et alJ (l2005h. that the PSF is constant 



or even, like lHevmans et al.1 (|2005f ). that the focus 
is piecewise constant. Fortunately, most of the 
PSF variations can be ascribed to a single physi- 
cal parameter: the distance between the primary 
and secondary mirrors or "effective focus" . Vari- 
ations of order 10/im create ellipticity variations 
of up to 5% at the edges of the field, which is 
overwhelming i n terms of a weak lensing signal. 
Jee et al.1 (|2005l ) built a PSF model for individual 



exposures by linearly interpolating between two 
PSF patterns, observed above and below nominal 
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focus. We have used the TinyTim (jKristl 120031 ) 
raytracing package to continuously model the PSF 
as a function of effective focus and CCD position. 
By matching the dozen or so stars brighter than 
FSUWab = 23 on eac h typical COSMOS image 
(|Leauthaud et all [2007) ;o TinyTim models, we 
can robustly estimate the offset fro m nominal fo- 



additive (c) shear errors. 

2.3. Charge transfer effects 



cus w ith an rms error less than 1/im ([Rhodes et al 



2007). We then return to the entire observational 
dataset, and fit a 3 x 2 x 2 order polynomial for 
each parameter of the PSF model, as a function of 
x, y, and focus. Using the entire COSMOS dataset 
strengthens the fit, especially at the extremes of 
used focus values, where few stars have been ob- 
served. The final PSF model for each exposure is 
then extracted from the 3D fit, at the appropriate 
focus value. 

We use the shear measurem ent method devel- 
oped f or space-based imaging bv lRhodes. Refregier k. 
(|2000L hereafter RRG). It is a "passive" method 
that measures the Gaussian-weighted second or- 
der moments ijj = ^^J}^ of each galaxy and 
corrects them using the Gaussian-weighted mo- 
ments of the PSF model. RRG is well suited to 
the small, diffraction-limited PSF obtained from 
space, because it corrects each moment individu- 
ally, and only divides them to form an ellipticity 
at the final stage. 

In an advance from previous implementations 
of KSB, and spurred by the fin dings of the 
Shea r TEsting Program (STEP; iMassev et al. 



20071 ). we allow the shear responsivity factor G 
to vary as a function of magnitude. The shear 
responsivity is the conversion factor between mea- 
sured galaxy ellipticity and the cosmologi- 
cal ly interesting quantity s hear 7.;. As described 



Leauthaud et alj (|2007l ). we have tested our 



pipel ine on simulated im ages created with the 
same lMassev et al.l ( 20041 ) package used for STEP, 
but tailored specifically to the image characteris- 
tics of the COSMOS data. We found it necessary 
to multiply our shears by a mean calibration factor 
of (0.86) -1 , but then found the shear calibration 
(m) accurate to 0.3%, with a residual shear offset 
(c) = 0.2 ±4 x I0 -4 , with no significant varia- 
tion as a function of simulated galaxy size or flux. 
This is particularly important in the measure- 
me nt of a shear signal a s a function of redshift . 
See lHevmans et all (|2006l ) or lMassev et~aH (|2007l ) 
for the definitions of the multiplicative (m) and 



As discussed further in iRhodes et al.l (|2007l ). 
the ACS WFC CCDs also suffer from imperfect 
charge transfer efficiency (CTE) during readout. 
This causes flux to be trailed behind objects, 
spuriously elongating them in a coherent direc- 
tion that mimics a lensing signal. Furthermore, 
since this effect is produced by a fixed number 
of charge traps in the silicon substrate, it affects 
faint sources (with a larger fraction of their flux 
being affected) more than bright ones. Thus it is 
an insidious effect that also mimics an increase in 
shear signal as a function of redshift. CTE trail- 
ing is a nonlinear transformation of the image, and 
prevents traditional tests of a weak lensing anal- 

vsis that look at bright stars. As such, it is the 

ehotft ost significant hurdle to overcome in weak lens- 
ing analysis from space. 

We are developing a method to remove CTE 
trailing at the pixel level. Following the work of 
Bristow ( 20041 ) on STIS, this will push charge back 



to where it belongs, as the very first stage in data 
reduction. Because an ACS version of this algo- 
rithm is still under development, in this paper we 
correct most of the CTE effect via a parametric 
model acting at the catalog level. We assume 
that the spurious change in an object's appar- 
ent ellipticity e, is an additive amount that de- 
pends only upon the object's flux, distance from 
the CCD readout register, and date of observa- 
tion. In fact, we also allowed variation with object 
size, although this h ad little effect. As shown in 



Rhodes et al. (|2007l) . this correction is sufficient 



for the full catalog of more than 70 galaxies per 
arcmin 2 when considering mass reconstruction, 
or circularly averaged statistics on small scales, 
where the signal is strong. However, it is not ad- 
equate for the faintest galaxies, when considering 
statistics on large scales, as we would like to do 
in this paper. Fortunately, the galaxy flux level 
at which the CTE correction successfully removes 
the CTE signal (leaving a residual signal one or- 
der of magnitude below the expected cosmological 
signal), appears to coincide with that for which 
reliable photometric redshifts can be obtained for 
almost all objects. 
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2.4. Photometric redshifts 

Reliable photometric redshift estimation is vital 
to the success of our 3D shear measurement. For 
this reason, the COSMOS field has been observed 
from the ground in a comprehensive range of wave- 
lengths (Capak et al. 2006). Deep imaging is cur- 
rently available in the Subaru Bj, Vj, g + , r + , i + , 
z+, NB816, CFHT it*, i*, CTIO/KPNO K s , and 
SDSS it', g', r\ i', z' bands. The COSMOS pho- 
tometric redshift code was used as described in 
Mobasher et al. (2006). This code contains a lu- 
minosity function prior in order to maximise the 
global accuracy of photometric redshifts for the 
faintest and most distant population. It returns 
both a best-fit redshift and a full redshift prob- 
ability distribution for each galaxy. The size of 
68% confidence limits for each estimated redshift 
are well-modelled by 0.03(l + z) o ut to z ~ 1.4 and 
down to magnitude Ifxuw = 24 (jMobasher et al 
2007t iLeauthaud et al.il2007l) . 



Before a large spectroscopic redshift sample be- 
comes available to calibrate the galaxy redshift 
distribution, our 3D analysis will be limited by 
the reliability of photometric redshifts. We do 
not impose a strict magnitude cut in the single 
Ifsuw band, but instead using color informa- 
tion from many bands, and select those galaxies 
with accurately measured redshifts. This includes 
96% of detected galaxies brighter than Ifsuw = 
24, and an incomplete sample fainter than that 



(jLeauthaud et al.l 120071 ) . The selection function, 
and the final redshift distribution, thus depend 
upon the spectral energy distribution of individual 
galaxies. However, since the background galaxies 
are unrelated to the foreground mass that is lens- 
ing them, such incompleteness has no detrimental 
effect on our analysis. 

We specifically select galaxies that are observed 
in the multi-color ground-based data and that 
have a 68% confidence limit in their redshift proba- 
bility distribution function smaller than Az = 0.5. 
The latter cut primarily removes galaxies with 
double peaks in the photometric redshift PDF due 
to redshift degeneracies. With the range of colors 
currently observed in the COSMOS field, one par- 
ticular degeneracy dominates: between 0.1 < z < 
0.3 and 1.5 < z < 3.2, where the 4000 A break can 
be confused with coronal line absorption features. 
At z > 1.5, the 4000 A break is well into the IR, 



where sufficiently deep data are not yet available 
for conclusive identification. To avoid catastrophic 
errors between these specific redshifts, we there- 
fore also exclude galaxies with any finite probabil- 
ity below z = 0.4 and above z = 1.0. After these 
cuts, we have redshift (and shear) measurements 
for 40 galaxies per arcmin 2 . 

3. 2D shear analysis 

3.1. 2D source redshift distribution 

The distribution of galaxies with reliably mea- 
sured shears and redshifts is shown in figure [1] 
The effects of cosmic variance are quite apparent, 
with the spikes below z ~ 1.2 all corresponding 
to known structures in the field. Beyond that, 
the photometric redshifts are limited by the fi- 
nite number of observed colors for each galaxy, 
and the peaks at z = 1.3, 1.5 and 2.2 arise arti- 
ficially at locations where spectral features move 
between filters. The median photometric redshift 
is Zmcd = 1-26. To minimize the impact of galaxy 
shape measurement noise, we downweight the con- 
tribution to the measured signal from faint and 
therefore noisier galaxies. We apply a weight 

W = —, \ , (1) 

<x £ (mag)+0.1 V ; 

where the rms dispersion of observed galaxy ellip- 
ticities is well-modelled by 



cr e (mag) w 0.32 + 0.0014(mag - 20) s 



(2) 



The error distribution of the shear estimators 



is discussed in more detail in Leauthaud et al 



(|2007t ). After this weighting, the median photo- 



metric redshift is 



1.11. In most cosmic 



shear analyses to date, an estimate of this value is 
all that was known about the redshift distribution. 
The smooth, dotted curve shows the distribution 



that w ould have been obtained from a lSmail et al 
( 19941 ) fitting function 



P(z) 



oc z 



exp 



(1.41z/z mod 



(3) 



with a = 2, j3 = 1.5, z me d = 1-26 and an over- 
all normalization to ensure the correct projected 
number density of galaxies. This would have been 
a better fit to the high redshift tail apparent in fig- 
ure [U had the free parameter in the model, z me( j, 
been ~ 1.17. 
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Fig. 1. — The thin, solid line shows the distribu- 
tion of the best-fit redshifts returned by the COS- 
MOS photometric redshift code (Mobasher et al. 
2006) with a luminosity function prior. The thick, 
solid line shows the distribution after accounting 
for the different weights given to galaxies. In both 
cases, the bin size is Az = 0.02. Peaks below 
z~ 1.2 correspond to real structures in the field, 
but the artificial clustering at higher redshift is 
due to limitations in the finite number of observed 
near-IR colors. The dashed curve shows the red- 
shift sensitivity function, assuming a ACDM uni- 
verse with WMAP parameters. The dotted line 
shows the redshift distribution that would have 
been expected, with knowl edge of only the me dian 
photometric redshift and a lSmail et al] (| 1994T ) fit- 
ting function. 

Figure [T] also shows the lensing sensitivity func- 
tion 



g(x) 



X \(r>) Da ^ Da ^' 
x V(X ' D A ( X ') 



^a-\ X ) d X ', 



(4) 

of the observed source redshift distribution, where 
X is a distance in comoving coordinates (in which 
the power spectrum is measured), \h is the dis- 
tance to the horizon, Das are angular diameter 
distances, (with the extra factor of a -1 converting 
these into comoving coordinates), and rj(x) is the 
distribution function of source galaxies in redshift 



space, normalized so that 

rXh 



v(x) d X = 1 



(5) 



o 



This represents the sensitivity of a projected lens- 
ing analysis to mass over densities, as a function of 
their redshift, and peaks at z ~ 0.4, about half- 
way to the peak of the source galaxy redshift dis- 
tribution in terms of angular diameter distance. 

3.2. 2D shear correlation functions 

The 2D power spectrum of the projected shear 
field is given by 



_9_ 

16 



Ho 

c 



Xh 



(6) 

where x is a comoving distance; Xh is the horizon 
distance; g(x) is the lensing weight function; and 
P(k,x) is the underlying 3D distribution of mass 
in the universe. The two-p oint shear correlations 
functions can be expressed (jSchneider et al.ll2002fl 
in terms of the projected power spectrum as 



cm 



4tt jo 
1 



CJ [3 o (e0) + 3 4 (£6)] ld£(7) 



C 2 (9) = — C] [J (£9) - J 4 (#)] I <W(8) 
4?r Jo 

These can be measured by averaging over 
galaxy pairs, as 

^{6) = <7i»7i> + 0) > (9) 
C 2 {9) = (l r 2 (r)Y 2 (r + e)) , (10) 

where 9 is the separation between the galaxies and 
the superscript r denotes components of shear ro- 
tated so that 71 (72) in each galaxy points along 
(at 45° from) the vector between the pair. In prac- 
tice, we compute this measurement in discrete bins 
of varying angular scale. However, they will need 
to be integrated later, so to keep this task man- 
ageable, we use fine bins of 0.1 "throughout the 
calculations, and only rebin for the sake of clarity 
in the final plots. 

A third shear-shear correlation function can be 
formed, 

C^{9) = < 7 [(r) l 2 \v + 9) ) + ( Y 2 (r) 7i>+0) ) , 

(11) 



G 



for which parity invariance of the universe requires 
a zero signal. The presence or absence of C^{9) can 
therefore be used as a first test for the presence of 
systematic errors in our measurement, although 
many systematics can still be imagined that would 
not show up in this test. 

The 2D shear correlation functions measured 
from the entire COSMOS survey are shown in 
figure [5J Note that the measurements on scales 
smaller than ~ 1' are new. For a given survey 
size, these are obtained more easily from space 
than from the ground because of the higher num- 
ber density of resolved galaxies. 

The additional, spurious signal that would have 
been obtained without correction for CTE trailing 
is shown as roughly horizontal, solid lines in fig- 
ure [21 This was calculated by recomputing the 
correlation functions but, rather than construct- 
ing a shear catalog by subtracting the CTE con- 
tamination from each galaxy's raw shear measure- 
ment, the CTE contamination was used as a di- 
rect replacement. An estimate of the residual CTE 
contamination for the galaxy population after cor- 
re ction, according to t he performance evaluation 



Rhodes et al.l (|2007l ). is shown as dotted lines. 



Although this is now below the signal, the uncor- 
rected level was more than an order of magnitude 
larger than the signal on large scales. Minimizing 
CTE by careful hardware design to avoid the need 
for this level of correction will be a vital aspect 
of dedicated space-based weak-lensing missions in 
the future. 

3.3. Error estimation and verification 

The error bars in figure [2] include both sta- 
tistical errors due to intrinsic galaxy shape noise 
within the survey, and the effect of sample ( "cos- 
mic") variance due to the finite survey size. The 
shape noise dominates on small angular scales, 
and the cosmic variance on scales larger than 
~ 10'. Surveys covering a similar area but in 
multiple lines of sight, s uch as ACS parallel data 
(jSchrabback et aill2007l Rhodes et al. in prepara- 
tion), will suffer less from the latter effect. 

The statistical shape noise is easy to measure 
from the gaalxy population. To measure the sam- 
ple variance, we split the COSMOS field into four 
equally-sized quadrants and recalculate the cor- 
relation functions in each. Of course, large-scale 




Fig. 2. — Correlation functions of the 2D shear 
field. The open circles indicate negative values. 
The inner error bars show statistical errors only; 
the outer error bars, visible only on large scales, 
also include the contribution of cosmic variance. 
The six parallel curves show theoretical predic- 
tions for a flat ACDM cosmology with fi TO = 0.3 
and erg varying from 0.7 (bottom) to 1.2 (top). 
The roughly horizontal lines indicate the level of 
the spurious signal due to CTE trailing before and 
after correction. 

correlations in the mass distribution mean that the 
four adjacent quadrants are not completely inde- 
pendent at large scales, and the measured variance 
underestimates the true error. To correct for this 
effect, we artificially increased the measured er- 
rors on 20-40' scales by 15%, in line with initial 
calculations. 

After the fact, we have compared our final er- 
ror bars to independent predictions from a full 
raytracing analysis through n-body simulations 
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Fig. 3. — Comparison of the error bars that we 
meas ured from the da t a, to advance predictions 
from Semboloni et al.l ( 2006h . obtained by ray- 
tracing through n-body simulations of large-scale 
structure. The two solid lines show the predic- 
tions assuming a Gaussianised mass distribution 
(bottom) and with the full, non-Gaussian distri- 
bution (top). 



by ISemboloni et al. ( 2006^ 1 . Figure [3] shows the 
predicted and observed la errors on C+(6) = 
Ci(8) + 6*2(6') (assuming 40 background galax- 
ies per square arcminute in the simulations, dis- 
tributed in redshift with z mc d = 1.11, and with 
<7 e = 0.32). Averaging across all thirteen angu- 
lar bins with equal weight, the mean ratio be- 
tween our measured error and the predicted non- 
Gaussian error is 0.994. Future work may there- 
fore improve the error estimation, but in the COS- 
MOS field at least, our quadrant technique reaches 
a level of precision sufficient for this paper. 

We also use the quadrant technique to measure 
the full covariance matrix between each angular 
bin. As shown in figure [4] the off-diagonal ele- 
ments are non-zero. This is expected even in an 
ideal case, because the same source population of 
galaxies is used to construct pairs separated by dif- 
ferent amounts. Nor are the upper-left and lower- 
right quadrants of figure[4]expected to be zero: the 
same pairs go into the calculation of both Ci(9) 
and C-2(9); and after deconvolution from the PSF, 



7[ and % are 110 longer formally independent. We 
will use the full, non-diagonal covariance matrix 
during our measurement of cosmological parame- 
ters in SJS] 

The final datum in the C${ff) panel of figure [2] 
is significantly (~ 5er) non-zero. This may be 
real: a finite region may not be parity invariant 
on scales comparable to the field size. But even if 
this does indicate a systematic problem, it is not 
as troubling as it appears, because on this scale 
the error bars are large for C\{9) and 62(6*), so 
the point carries very little weight. For a possi- 
ble explanation, note that the spurious C3 signal 
has the same sign as the uncorrected CTE signal. 
On scales that span almost the entire COSMOS 
survey, one of the galaxies in a pair must lie near 
the edge of the survey field - which was observed 
last, and suffers most from CTE degradation. If 
the temporal dependence of the CTE signal is not 
linear, as we have assumed, the spiral observing 
strategy to cover the field would produce a simi- 
lar CTE pattern (and a coherent residual signal) 
in all four quadrants. This could create an ad- 
ditional 63(40') signal, with error bars underes- 
timated by our quadrant method. Resolving this 
issue requires CTE data from a longer time span, 
or more data separated by 20-40'. Such analysis 
may be feasible with ACS parallel data (Rhodes 
et al. in preparation) but is not possible here. 

3.4. 2D shear variance 

For historical reasons, cosmic shear results are 
often expressed as the variance of the shear field 
in circular cells on the sky. For a top-hat cell of 
radius 9, this measure is related to the shear cor- 
relation functions by 
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dtf, (12) 



where we have used a small angle approximation. 
Note that the signal is more strongly correlated on 
different angular scales in this form than it is when 
expressed as correlation functions. The results are 
shown in figure [5] 

3.5. 2D E-B decomposition 

The correlation functions can also be recast in 
terms of non-local E (gra dient) and B (cur l ) pat- 
terns in the shear field ( Crittenden et al. 2000; 
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Fig. 4. — Covariance matrix for the 2D correla- 
tion functions C\ (9) and C2 (9) shown in figure [2j 
obtained by splitting the COSMOS field into four 
quadrants and performing the analysis separately 
in each. The diagonal elements illustrate the size 
of the errors in each of the thirteen 9 bins, and 
the off-diagonal elements illustrate how much the 
measurements are correlated. The color scale is 
logarithmic. 



Pen et al]|2002l ). Gravitational lensing is expected 
to produce only E modes, except for a very low 
level of B mod es due to lens-lens co upling along 
a line of sight (jSchneider et al.ll2002l ). It is com- 



monly assumed that systematic effects would af- 
fect both E- and _B-modes equally. The presence 
of a non-zero B mode is therefore a useful indica- 
tion of contamination from other sources. 

E- and i?-modes correspond to patterns within 
an extended region on the sky, and cannot be 
separated locally. As a result, this operation for- 
mally requires an integration of the shear correla- 
tion functions over a wide range of angular scales. 
Two mathematica l funct i ons have been develope d 



((Crittenden et all l200Cj ISchneider et~aT1 [2002), 



which each include an in tegral over only small 
scales or large scales (see ISchneider fc Kilbingerl 
( 20061 ) for a new suggestion to construct a third). 
However, neither integral is ideal in practice, be- 
cause our correlation functions are only well mea- 
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Fig. 5. — Variance of the 2D shear signal in cir- 
cular cells of varying size. Solid lines show predic- 
tions in a concordance cosmology with as varying 
as in figure [2] Note that adjacent data points are 
highly correlated. 

sured on scales between ~ 0.5' and 40'. The 
absence of complete data introduces an unknown 
constant of integration, and it is not possible to 
uniquely split this measured shear field into dis- 
tinct E- and B-mode components. As a practical 
attempt to estimate this constant, we extrapolate 
data into the unknown regime, using predictions 
from the best-fit cosmology that is determined in 

m 

The signal on large angular scales is small, and 
the corresponding integrals require the least cor- 
rection. To calculate these, we first define C+ = 
C\ + C2 and C_ = C\ — C2. Then we can compute 



i E {Q) = C x {9) 



1 - 



39 2 \ C_(tf) 







dtf , 
(13) 



which contains only the E-mo&e signal and 



i B {9) = C 2 (9) 



1 - 



36» 2 \ C-(d) 



2 







di) 



(14) 

which contains only the S-mode signal. It is gen- 
erally necessary to add a function of 9 (not only 
a constant of in tegration) to £.f,( 9) and subtract it 
from £ B (6) (c.f. iPen et al~ll200^ . 
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Fig. 6. — E-B decompositions of the 2D shear 
field. The top panel shows the statistics that for- 
mally require an integral over the measured corre- 
lation functions to infinite scales, and the bottom 
panel shows those that formally require an integral 
from zero. Solid points show the -E-mode, and 
open points show the B-mode. The lines show 
predictions in a concordance cosmology with as 
varying as in figure [2 Note that adjacent data 
points arc in the top panel are highly correlated. 

The components can also be separated via the 



variance of the aperture mass statistic M ap (9). 
This is obtained from a weighted mean of the tan- 
gential (74) and radial (7,.) components of shear 
relative to the center of a circular aperture. This 
statistic is given by 

M &P {6) = J w(\$\;t)) 7t (tf) d 2 tf (15) 

which contains only contributions from the E- 
mode signal and 



M ± (0) 



W 



(\d\-e) lr (d) d 2 tf, (16) 



which contains only the B-mode signal, where 
W(|i9|) is a compensated filter. We adopt a com- 
pensated "Mexican hat" weight function 



W (^0) = ^^(l--)lht> ■,)). (IT) 



2 



where 8 defines an angular scale of the aperture 
and the Heaviside step function H truncates the 
weight function on large scales. 



Schneider et al.1 (|2002l ) derived expressions for 



the variance of these statistics as the aperture is 
moved across the sky. These reqire integrals over 
the correlation functions from small scales 



26 d_ 

(9 2 



(Ml) (0) 



211 



2 



where 



T+(x) = 



6(2- 15x 2 ) 



C_(tf)T_ 

C+WT+ 
C_(#)T_ 



1 arcsin 

7T 



dtf (18) 



dtf, (19) 



(!) 



(20) 
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100tt 
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Fig. 7. — The thin, solid line shows the redshift 
distribution of source galaxies. The thick, solid 
line shows their distribution after accounting for 
the magnitude-dependent weighting scheme. In 
both cases, the bin size is Az — 0.02. The dashed 
lines show (artificially normalized) redshift sensi- 
tivity curves obtained by slicing this distribution 
into the discrete redshift bins indicated by the ar- 
rows at the top. 

for x < 2 and T+(x) = T-(x) = for x > 2. We 
again estimate the constant of integration by ex- 
trapolating our data with theoretical predictions 
in cosmological model preferred by the rest of the 
data. 

From figure [H we can see that £b(9) is consis- 
tent with zero on all scales. The noise is particu- 
larly large on small scales, and the rather unstable 
M\_(ff) is affected on scales up to ~ 1' by the first 
bin. 

4. 3D shear analysis 

4.1. Correlation function tomography 

We now split the catalog into three discrete red- 
shift bins and, as before, calculate the correlation 
functions using all pairs of galaxies within each 
bin. The redshift bins are chosen in considera- 
tion of the particular color information available. 
Degeneracies in the photometric redshift estima- 
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Fig. 8. — Evolution of the cosmic shear two-point 
correlation function signal with increasing red- 
shift. The series of data points, from bottom to 
top, show measurements from slices between red- 
shifts 0.1, 1, 1.4, and 3. The black curves show pre- 
dictions from a flat ACDM model with f2 m = 0.3 
and erg = 0.85, for the same slices, increasing in 
redshift from the bottom to the top. Open circles 
depict negative values. 

tion cause galaxies with a flat distribution in red- 
shift to cluster artificially around z = 1.3, 1.6 and 
2.2. An excess at these positions is evident in 
figure [T] We therefore pick bins with boundaries 
away from these values, and with widths similar to 
the size of the local peaks in the resdshift distribu- 
tion. For COSMOS, suitable bins are 0.1 s$ z sC 1, 
1 < z ^ 1.4 and 1.4 < z ^ 3. This scheme conve- 
niently divides up the galaxies fairly evenly, with 
the slices each containing 32%, 24% and 44% of 
the galaxies. Unfortunately, the last bin can not 
be further sub-divided without deeper IR or UV 
data. The redshift slices and their resulting lens- 
ing sensitivity functions are illustrated in figure [7] 
Figure [8] shows the increasing two point cor- 
relation function signal for pairs of source galax- 
ies as a function of redshift, where both galaxies 
are in the same redshift bin. Since the measure- 
ments in the redshift bins are much more noisy 
than those from the projected 2D analysis, we plot 
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C+(9) = C 1 {0)+C 2 (0) infigure[8] Theoretical pre- 
dictions for the correlation functions are obtained 
for each slice by replacing the lensing weight func- 
tion g(z) in equation © by those shown in fig- 
ure and obtained from only the galaxies in a 
given slice. Because the effective lensing volume 
J g(z) dz increases for succesive redshift bins, the 
signal increases with z. 

Figure [5] shows the measured covariance ma- 
trix for the 3D correlation functions. The de- 
gree of correlation between the lowest and highest 
redshift bins, primarily evident on small scales, 
is unexpected. Had it been significant on all 
scales, a likely explanation would have been cross- 
contamination of the bins by galaxies from other 
redshifts (the well-known degeneracy between low 
and high redshift from photo- z estimation is dis- 
cussed in £ )2.4[) . Had the covariance been equally 
evident in all three bins, likely explanations 
could have been: interference o f intrinsic align- 
ment s like those suggested by (jHirata fc Seliakl 
2004); and imperfect correction for PSF varia- 
tion or DRIZZLE-related pixellisation effects un- 
accounted for on small scales. In practice, the 
most likely explanation is a combination of sev- 
eral such effects, each at a low level. 

Although the signal in the individual slices is 
noisy, we have attempted an E-B decomposition 
in figure [TOJ using the same two statistics as were 
applied to the 2D analysis. The integrals over 
the noisy correlation functions are particularly ill- 
defined at < V , Nevertheless, the signal in- 
creases to high redshift, matching the theoretical 
expectation for this measurement. 

4.2. Growth of structure 

The total i?-mode signal corresponds to the 
integrated mass density along a line of sight, 
weighted by the lensing sensitivity function. The 
evolving i?-mode signal in figures [5] and [TO] grows 
towards high redshift due to the increasing volume 
that it probes, and in which mass structures are 
located. This offers constraints on the large-scale 
geometry of the universe. But, if we are more 
interested in the mass structures themselves, this 
function of in fixed redshift slices can be recast 
into a function of z for fixed angular scales. We 
shall now suggest a new way of viewing this data, 
which stays close to measurable quantities, but 
offers a new insight into the underlying structure 
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Fig. 9. — Covariance matrix of correlation func- 
tion data used in the 3D cosmic shear analysis. 
Note that this includes only the three correlation 
functions where both galaxies are in the same red- 
shift bin. An additional three correlation functions 
can be formed from pairs in which the galaxies 
come from different slices, but these are not shown 
in this plot for the sake of clarity. 

formation. 

Each data point in figure [8] corresponds to the 
amount of mass within an effective volume. This 
volume is described in azimuthal directions by 
Bessel functions, and in the redshift direction by 
the lensing sensitivity function g(z). Assuming 
the best-fit cosmology from J}5] to fix the geome- 
try of the universe, we can divide by this volume, 
and obtain a quantity proportional to mass den- 
sity. In practice, to increase the signal-to-noise of a 
measurement that will involve many redshift bins, 
we do not restrict the measurement to only those 
pairs within a given redshift slice, as before. We 
require the nearer galaxy to be inside the slice, 
but then compute correlation functions using all 
galaxies behind it. The more distant galaxy has 
then been lensed by anything the foreground has 
been lensed by. The effect is merely to change the 
(squared) lensing sensitivity function to the prod- 
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Fig. 10. — E-B decomposition of the 3D cosmic 
shear signal, in different redshift bins, colored as in 
figure [51 For clarity, only the -E-modes are shown. 
Open circles depict negative values. The f?-modes 
are as noisy, but are consistent with zero. Note 
that adjacent data points are highly correlated. 



uct of the sensitivity function for the slice galaxies 
with that of the background distribution. This 
creates a new, effective g(z) that peaks at slightly 
higher redshift, but is still zero behind the nearest 
galaxy. 




Age of the universe [Gyr] 



Fig. 11. — The growth of structure over cosmic 
time. This links the cosmic shear signal on fixed 
angular scales as a function of redshift (rather than 
the other way round, as in previous figures). Data 
points are located at the peak of the lensing sen- 
sitivity function for each set of source galaxies. 
The source galaxies themselves are approximately 
twice as far away. The different colors distinguish 
different angular scales. For each of these, the 
dashed line shows the theoretical expectation, as- 
suming the best-fit cosmological model from <J5] 



Figure QT] thus shows 



G(z;8) = C+(6;z) / f g 2 {z') dz' (22) 



1 JC](d;z) J (£6) m 
2^ Jq g 2 (z') dz' 



(23) 



the growth of power on different angular scales. 
The foreground mass is most likely to lie near the 
peak of the sensitivity function, so we place the 
data points at this redshift. In practice, it could 
lie anywhere within g(z), so we overlay error bars 
in z equal to the rms of g about its peak. The- 
oretical predictions of this quantity are overlaid, 
assuming a flat, ACDM cosmology, with the best- 
fit parameters found in Sj^l 

The growth towards z = represents a com- 
bination of the physical growth of structure and 
the mixing of fixed physical scales at different red- 
shifts into a measurement at one apparent angu- 
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lar scale. Both of these effects act in the same 
sense, to increase the signal towards the present 
day. This is in contrast to the cosmic shear signal 
in figures [5] and [TUl which itself increases towards 
high rcdshift. On large scales, the small cosmic 
shear signal makes the measurement fairly noisy. 
On intermediate scales, the data closely follow the 
predictions. The lowest redshift point is obtained 
from pairs of galaxies where the nearest is between 
z = 0.1 and z — 0.7. We speculate that the ap- 
parently significant upturn at low z and on small 
scales might be caused by contamination of that 
redshift bin by high redshift galaxies. These could 
have been caught by the photometric redshift de- 
generacy discussed in H'2A\ and would contain a 
apparently spurious signal when moved to low red- 
shift. The accuracy of the photometric redshifts 
may therefore be limiting the precision of this mea- 
surement. 

5. Constraints on cosmological parame- 
ters 

5.1. 2D parameter constraints 

We now use a Maximum Likelihood method to 
determine the constraints set by our 2D obser- 
vations of Ci(9) and C^f) upon the cosmologi- 
cal parameters f2 m , the total mass-density of the 
universe, and erg, the normalization of the matter 
power spectrum at 8 hr 1 Mpc. We assume a fiat 
universe, with a Hubble parameter h = 0.7. 



We closely follow the approach of jMassev et al 



(|2005l ). obtaining theoretical predictions for the 
lin ear transfer function f rom the fitting functions 



of Bardecn et al 



power spect r um 



Smith et al 



1 1986f ) and for the non- linear 
using the fitting functions of 



( 20031) . The theoretical correlation 
functions are first calculated from equation ([6]) 
in a three dimensional grid spanning variations 
in fl m from 0.05 to 1.1, as from 0.35 to 1.4 and 
the power spectrum shape parameter T from 0.13 
to 0.33. We used the full redshift distribution 
of source galaxies (after correction for weighting) 
shown in figured] 

We then fitted the observed shear correlation 
functions d(9) — {Ci(9), 6*2(6')} to the theoretical 
predictions calculated at the centers of each bin 
t(9), computing the log- likelihood function 

X 2 (tf, Q m , a s ,T) = (d{9) - t{9, Q m , a 8 ,T)) T cov(d)' 




Fig. 12. — Constraints on cosmological parame- 
ters from a traditional 2D cosmic shear analysis, 
after marginalization over other free parameters. 
In order of decreasing thickness, the contours in- 
dicate 68.3%, 95.4% and 99.7% confidence lim- 
its due to statistical errors; additional uncertainty 
potentially contributed by sources of systematic 
error are discussed in the text. The grey scale 
background is logarithmic, and shows y 2 divided 
by the number of degrees of freedom in the data. 
The white a rea at the bo t tom-r ight was excluded 
because the ISmith et al.l ([2003) fitting functions 
could not be evaluated without unreasonable ex- 
trapolation of the non-linear matter power spec- 
trum to physical scales smaller than 0.1/i _1 kpc. 
This can be compared to the much tighter con- 
straints from the full 3D analysis in figure 1151 

throughout the grid, where cov(d) is the covari- 
ance matrix in figure [¥] In an advance of earlier 
incarnations, we perform the matrix inversion via 
a singular value decomposition (SVD), and discard 
all eigenvalues not within machine precision of 
the largest. We do n ot include the multiplicative 
factor suggested by lHartlap. Simon &: Schneider 
([20071 ) . We then marginalize over T with a Gaus- 
sian pri or centered on 0.19 an d with an rms width 
of 15% ( Percival et al.ll200lh . To compute confi- 
dence contours, we numerically integrate the like- 
lihood function 



1 (d(0)-tt0,fi m ,<T 8 ,r)) 



-X 2 /2 



(24) 
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Our constraints on cosmological parameters 
from this 2D analysis are presented as a pro- 
jection through parameter space in figure 1121 
The contours represent statistical errors, includ- 
ing full non-Gaussian sample variance. Formally, 
the best-fit model has O m = 0.30, a% — 0.81 
and T = 0.21, and this achieves Xreduccd = 
X 2 /( n param — 3) = 1.10 in 23 degrees of free- 
dom. However, there is a well-known degeneracy 
between f2 m and as when using only two-point 
statistics. Changing the T parameter slides the 
contours back and forth along this valley, and 
marginalization over this parameter also slightly 
increases the minimum \ 2 - After marginaliza- 
tion, a good fit to our 68.3% confidence level from 
statistical errors is given by 
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"03 



0.44 



0.81 ±0.075 



(25) 



with 0.15 fl m ^ 0.7. 



Massev et al.l (|2005j ) were unable to use the full 



covariance matrix due to instabilities in the ma- 
trix inversion, so had set to zero any elements in 
the covariance of C\ (6) with Ci (9) (these are the 
bottom- left and top-right quarters in figure [4j. 
This problem has been resolved in the present 
work by the use of an SVD. However if we discard 



half o f the covariance matrix as in iMassey et al 



(120051) . we obtain parameter constraints 



^8 



03 



0.44 



= 0.83 ±0.07 



(26) 



If we discard all of the off-diagonal elements in the 
covariance matrix, we obtain 



0"8 



03 



0.84 ±0.065 



(27) 



The slightly smaller error bars are expected, but 
the shift in the best-fit value relative to result ([25]) 
is not. This effect might go some way towards 
explaining the hig her than usua l value obtained 
for this quantity in iMassev et al. ( 2005 ). 

Note that all of the above constraints incor- 
porate only statistical sources of error; although 
these do include non-Gaussian sample variance 
and marginalization over other parameters. We 
can propagate the various sources of potential sys- 
tematic error by noting that 



C;(5') CX n 



1.46 Jl. 45 „1.65t 



-0.11 



(P 7)- 



(28) 



for i S {1,2} in a fiducial ACDM cosmological 
model with Sl m = 0.3, S1 A =0.7, T = 0.21 and 
o"8 = 1-0. Adding an uncertainty equivalent to 
10% in the median source redshift, a 6% shear cal- 



ibration uncertainty ( see [L cautha ud et alj 12007 



Hevmans et al.l 120061 IMassev et al.l 120071) . and 



an empirically esti mated binning instability (c.f. 
Massey et al.ll2005l ) to our constraint from the full 
covariance matrix gives a final 68.3% confidence 
limit of 



08 



03 



0.48 



: 0.81 ± 0.075 ± 0.024 ± 0.05 ± 0.02 
= 0.81 ±0.17, (29) 



where the various systematic errors have been 
combined linearly on the second line. 

5.2. 3D parameter constraints 

We shall now include the redshift information 
available for each object, adopting the 3D binning 
scheme introduced in <gl A simple 2D analysis 
can first be performed within each redshift slice, 
by simply exchanging the redshift sensitivity func- 
tion g(z) calculated using the full redshift distri- 
bution for one calculated using the restricted dis- 
tributions. Figure [TBI shows the constraints on cos- 
mological parameters from each slice, using only 
pairs of galaxies where both pairs lie in that slice, 
but the full covariance matrix for each. The in- 
dividual results are clearly more noisy than for 
the full 2D analysis, since each slice contains only 
approximately one ninth of the number of galaxy 
pairs. However, all of the slices are consistent with 
our base cosmological model. Furthermore, while 
the statistical noise is similar in each slice, because 
they all contain a similar number of galaxy pairs, 
the signal (and hence the signal to noise) clearly 
increases at high redshift, as expected. 

In figure [JJJ the constraints from the three red- 
shift bins are combined as if they all provided inde- 
pendent information (despite the fact that the red- 
shift sensitivity functions in figure [7] clearly over- 
lap, so they are correlated). Although there are 
approximately only one third of the number of 
galaxy pairs in this analysis as there were in the 2D 
analysis, the additional information about the evo- 
lution of the signal as a function of redshift retight- 
ens the 68% confidence limit constraints back to a 
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Fig. 13. — Constraints on cosmological parame- 
ters from within each of the three separate redshift 
slices, from low (top panel) to high redshift (bot- 
tom panel) . The redshift binning scheme is shown 
in figure [7] and discussed in the text. The contours 
indicate 68.3%, 95.4% and 99.7% confidence lim- 
its, and the logarithmic color scale is common to 
all three slices. 



10' 



10* 




Fig. 14. — Combined constraints on cosmological 
parameters f2 m and as from a series of effectively 
2D shear analyses in each of the three redshift 
slices (see text). Only pairs of galaxies where both 
lie in the same redshift slice have been included in 
this analysis. This can be compared to the similar 
result from the 2D analysis in figure Q21 and the 
full 3D analysis in figure [j"5l 



similar value of 
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03 



0.44 



= 0.86 ±0.08 



(30) 



for Q m ^ 0.25. The best-fit model has Q m = 0.55 
and c 8 = 0.64, which achieves X reduccd = 1-18 in 
28 degrees of freedom. 

We can restore the missing galaxy pairs, and 
their information content, by introducing three 
additional correlation functions, constructed from 
pairs of galaxies that lie in different redshift slices. 
The theoretical expectation for these correlation 
functions requires the g 2 (z) term in equation ([6]) 
to be replaced by the product of the lensing sen- 
sitivity functions for the two redshift bins. We 
use the full covariance matrix, which is again esti- 
mated from variation between the four quadrants 
of the COSMOS field. Figure [15] shows a projec- 
tion of the log-likelihood surface, with the usual 
contours. 

The best-fit model has ft m — 0.47 and ug = 
0.72, which achieves Xreduced = 2.35 in 56 de- 
grees of freedom. This is significantly greater than 
unity because only statistical errors are currently 
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Fig. 15. — Constraints on cosmological parame- 
ters Q m and o~%, from a full 3D cosmic shear anal- 
ysis. Solid contours indicate 68.3%, 95.4% and 
99.7% confidence limits due to statistical errors 
and marginalization over other parameters; po- 
tential sources of additional, systematic error are 
discussed in the text. These constraints are far 
tighter than the equivalent results from our simple 
2D analysis, which are reproduced from figure [T2l 
as dotted lines for ease of comparison. The white 
area at the bottom -right was excluded because the 
Smith et all {2003) fitting functions could not be 
evaluated without unreasonable extrapolation of 
the non-linear matter power spectrum to physical 
scales smaller than 0.1/i _1 kpc. 



included. As described below, the error budget is 
increased by a factor of ~ 1.5, and the minimum 
X?oducod to 1.04, when considering systematic er- 
rors in the relative shear calibration and mixing 
of galaxies between bins. Again we find the usual 
degeneracy, the best-fit position along which is de- 
termined by the parameter T. However, with the 
full 3D information, parameter constraints in the 
direction orthogonal to this are much tighter. Our 
68% confidence limits are well-fit by 



(13 



0.44 



°8 ( tttt ) = 0.866 ± 0.033 , 



(31) 



for 0.3 s: Sl m < 0.6. 

We now incorporate a systematic error bud- 
get into our 3D parameter constraints. We al- 



low a 6% absolute she ar calibration uncertainty 
(jLeauthaud et al.l 120071 ). a 5% relative shear cal- 
ibration uncertainty between low and high red- 
shift bins, and a pot ential 10% contamination (e.g. 
iMassev et al. 2004 ) of the high redshift bin by 
galaxies really at low redshift (and vice versa) due 
to the possibility of catastrophic redshift errors 
discussed in £12.41 This leaves a final 68.3% confi- 
dence limit of 



£7 m 
03 



: 0.866 ± 0.033 ± 0.026 ± 0.009±g;^ 



= 0.866™, 



(32) 



where the various systematic errors have been 
combined linearly on the second line. Note that, 
when considering the relative improvement in the 
parameter constraints from a 2D analysis (|29| to 
a 3D analysis (|32[) . it is not appropriate to in- 
clude errors from uncertainty in the absolute cal- 
ibration of a shear measurement method that is 
common to both. Continuing to budget for po- 
tential relative mis-calibration between low- and 
high-redshift bins, as well as including all other 
sources of systematic and statistical error, reveals 
a dramatic threefold tightening of parameter con- 
straints. 

We have also tried increasing the number of red- 
shift slices, for a finer quantitative measurement of 
the evolution of the shear signal. We attempted an 
analysis using five redshift bins, created by split- 
ting in half the first two slices of the three used 
previously. Unfortunately, the covariance matrix 
became degenerate, and harder to invert. Fur- 
thermore, the best-fitXreduced an d cosmological pa- 
rameter constraints degraded. The results in each 
bin were very noisy (the signal to noise is pro- 
portional to ™g a 2 laxics ), but, as in H4.21 there were 
hints that the signal did not evolve as expected 
after this finer redshift binning. The likelihood 
surfaces from individual slices did not agree, so 
their combination was blurred out. We interpret 
this as indicating that galaxies were beginning to 
be placed in the wrong redshift bins, and polluting 
that signal. Thus we have effectively reached the 
available precision of the photometric redshifts, at 
least at the high redshifts in which the weak lens- 
ing signal is concentrated. For further progress, 
we await ongoing, deeper multicolor observations 
of the COSMOS field. 
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6. Conclusions 

We have performed a fully three dimensional 
cosmic shear analysis of the largest ever survey 
with the Hubble Space Telescope. The 3D shear 
field contains rich information about the growth 
of structure and the expansion history of the uni- 
verse. Indeed, by assuming a concordance cos- 
mological model, we have directly measured the 
growth of structure on both linear and non-linear 
physical scales. We have also placed independent 
68% confidence limits on cosmological parameters. 
From a traditional, two dimensional cosmic shear 
analysis, we measure 



<78 



£7 m 
03 



0.48 



= 0.81 ±0.17 



(33) 



with 0.15 ^ 17 m ^ 0.7. From a full, three dimen- 
sional analysis of the same data, we obtain 



CT8 



03 



0.866 



-0.085 
-0.068 



(34) 



with Q m ^ 0.3. This represents a dramatic im- 
provement over already remarkable constraints. In 
fact, disregarding uncertainty in the absolute cal- 
ibration of our shear measurement method, which 
is common to both analyses, the 3D constraints 
represent a threefold relative improvement in the 
errors from the 2D constraints. 

A solely two-point cosmic shear analysis cannot 
easily isolate a measurement of just fi m . The de- 
generacy with as is broken only by the difference in 
signal between large and small scales. Our best-fit 
value of £l m is slightly larger than the measure- 
ment of 0.23 ± 0.02 from the 2dF galaxy redshift 
survey (jCole et al.l 120051 ). As discussed in 



undetected CTE correction residuals could poten- 
tially affect our measurement on the very largest 
scale in each redshift slice; however this datum 
carries very little weight because of shot noise, so 
this explanation is unlikely. Because high redshift 
slices probe larger physical scales than low redshift 
slices, our measurement of f2 m could also be po- 
tentially biased by photometric redshift failures. 
After allowing for this effect in our systematic er- 
ror budget, the discrepancy in f2 m is within la of 
the 2dF results, so we shall not pursue this further. 

The main constraint from our data is upon 
as- We find a value slightly larger than that of 



0.74+ °^ from the 3-year WMAP data (jSpergel et al 
2006). Our result is also larger than most esti- 
mate s of cluster abundance from x-ray survey s 
(e.g. iBorgani et all |200ll ISchneider etal] 120021 ), 



and from othe r recent space-based weak lensing 



meas urements (jHevmans et al.l2 005; Sc hrabback et al 
20071 ). However, the HST GEMS survey, on 
which both of the latter were based, suffers 
from sample variance due to its limited size, 
and is suspected from other measures of con- 
taining an unusually empty portion of the uni- 
verse. Furthermore, independent measurements 
of as — 0.85 or sli ghtly greater have recently 
been published by iMcCarthv. Bower &; Baloghl 
(l2006l) . : rom observations of the gas mass frac 



tion in x-ray selected clusters; Li et al. ( 20061 ), by 

counting the number of observ ed g iant arcs; and 

Viel, Haenelt k Springel (|2004l ) and lSeliak. Slosar fc McDonald! 
( 20061 ), with Lvman-q forest data. All of these 



measures contain information about small-scale 
density fluctuations at relatively low redshift: 
something much more intrinsically suited to a 
measurement of as than the CMB. Our results 
are also remarkably consistent with those from the 
ground-based CFHT wide synoptic legacy survey 
( Hoekstra et al.l l2006f ) . Such agreement between 
the largest space-based and ground-based surveys 
demonstrates the maturity of the field post-STEP. 
The combination of all these results is therefore 
beginning to hint at inconsistencies in cither the 
standard cosmological model or in the interpreta- 
tion of one or more of these methods. 

With the profundity of this statement in mind, 
we are careful to realistically include all possible 
sources of systematic error. The dominant contri- 
bution to the total error budget is uncertainty in 
the absolute calibration of our shear measurement 
method. The weak lensing community is earnestly 
working to improve and ascertain the reliability of 
various methods through simul ated images that 
contain a known input signal dLeauthaud et all 
2007t iHevmans et al.ll200& iMassev et al.ll2007h 

Aside from this contribution, further exploita- 
tion of the COSMOS survey is currently limited by 
two additional sources of potential systematic er- 
ror. Conveniently, these two limits currently hap- 
pen to lie at a similar flux level and therefore affect 
a similar population of galaxies, which we simply 
remove from our analysis. Since a weak lensing 
measurement is concerned with the mass distri- 
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bution in front of galaxies rather than the galax- 
ies themselves, this can be done without worries 
about bias. Firstly, the in-orbit degradation of the 
ACS CCDs has led to inadequate charge transfer 
efficiency during readout, which creates trailing of 
fai nt objects, and mim ics a weak lensing signal. 
In lRhodes et al.1 we formulated an empir- 

ical correction scheme for the CTE effect, which 
works for all but the faintest galaxies; an ongoing 
effort to correct CTE pixel-by-pixel in raw images 
should allow us to push this limit and dramatically 
increase the number density of galaxies with mea- 
sured shears. Secondly, the finite number of col- 
ors available for each galaxy, and particularly the 
depth in near-IR bands, limits the current accu- 
racy of photometric redshifts. Continuing obser- 
vations with the Subaru telescope should improve 
their precision. This will allow finer resolution in 
the redshift direction and, most importantly, will 
break redshift degeneracies ubiquitous in the red- 
shifts of faint objects so that they can also be used. 

By understanding the characteristics of effects 
that dominate real data, COSMOS is proving an 
invaluable dry run for future, dedicated weak lens- 
ing missions in space. We have revealed important 
aspects that should ideally be minimized by hard- 
ware design and mission scheduling requirements. 
However, we have also demonstrated the rich in- 
formation content of the 3D shear field, and shown 
a proof of concept for some of the proposed tomo- 
graphic analysis techniques that will be required 
to fully exploit such future data. 
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